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In their comment to our Letter on aging in a Lennard- 
Jones glass Q] , Mussel and Rieger || give evidence that a 
different type of scaling function, specified below, might 
give rise to a better scaling of the two-time correlation 
functions C q (t w ,t w +t) than the scaling function propsed 
by us. By using data stemming from their own simula- 
tion, they show that a scaling of the form 

C q (t w ,t w + 1) ~ <7{ln((i + t w )/r)/ ln(^/r)}, (1) 

where r is a fit parameter, leads to a better collapse of 
the curves for the different waiting times t w , than the 
scaling used by us, which is of the form C q (t w ,t w + 1) ~ 
C{t/t r }, with t r oc t^, a « 0.9, and which thus almost 
corresponds to "simple aging" (t r oc t w ). We therefore 
used the scaling function proposed by Mussel and Rieger 
in order to see whether this function can be used to scale 
also the data from our simulation onto a master function 
and show the result in Fig 1. The value of r is 0.005, 
which is comparable to the one found by Mussel and 
Rieger for their data. As is evident from the figure, this 
sort of scaling (main figure) does not work well for our 
data and is clearly inferior than the scaling proposed by 
us (inset). Thus we conclude that the scaling function 
proposed by Mussel and Rieger is not always appropriate 
to scale the two-time correlation functions for this sort 
of system. 

What remains unclear for the moment is the reason 
why the scaling function given by Eq. (1) works for the 
data of Mussel and Rieger, whereas it fails to do so for 
our data. A careful comparison of their and our data 
shows that the two sets of curves show some systematic 
differences, in that, e.g, the height of the plateau at in- 
termediate times is slightly larger in their data than in 
our data. Also, as stated by Mussel and Rieger, the t w 
dependence of the relaxation time t r differs from ours, 
since they find that the exponent a is 1.1 (as oposed to 
0.88). The reason for the difference of the relaxation data 
might be that the two simulations were not carried out 
in exactly the same way. For example the size of the 
system is different (32768 vs. 1000), the time step is dif- 
ferent (0.01 vs. 0.02) and the value of the wave-vector 
is different (7.5 vs. 7.25). It is not clear which ones of 
these differences, if any, gives rise to the slightly different 
aging behavior. It has been found before, however, that 


the details of the nonequilibrium dyamics can depend 
quite sensitively on the details of the simulations [E10. 
Finally we mention that if Tf, the final temperature of 
the quench, is decreased significantly, the scaling behav- 
ior found in Ref. Q does not seem to hold anymore and 
that also the Ansatz given by Eq. (1) does not seem to 
work 101 


[1] W. Kob and J.-L. Barrat, Phys. Rev. Lett. 78, 4581 
(1997). 

[2] U. Mussel and H. Rieger, Phys. Rev. Lett, xx, xxxx 
(199x). 


[3] G. Parisi, Report No. |cond-mat/9701015 
[4] J.-L. Barrat and W. Kob (unpublished). 


1.0 


+ 


o 


t 0.8 
0.6 
0.4 
0.2 



0.0 


1 V' 

10 t/t jo 


-7 -fi Ls —d -1 -? -1 

10 10 10 10 10 10 10 


10 L 


ln((t+0/T)/ln(VT)-1 


FIG. 1. Main figure: The data of Ref. [1] scaled in the way 
proposed by Mussel and Rieger with r = 0.005. Note that 
the data for t m = is not shown. Inset: The same data but 
scaled by t T oc t°- 88 . 
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